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Our team accomplished the following results in the research areas: 


Godunov-type numerical solution of the 20 moment plasma transport equa- 
tions. One of the centerpieces of our proposal was the development of a 
higher order Godunov-type numerical scheme to solve the gyration domi- 
nated 20 moment transport equations. In the first step we explored some 
fundamental analytic properties of the 20 moment transport equations for 
a low b plasma, including the eigenvectors and eigenvalues of propagating 
disturbances. The eigenvalues correspond to wave speeds, while the eigen- 
vectors characterize the transported physical quantities. In this paper we 
also explored the physically meaningful parameter range of the normalized 
heat flow components. In the second step a new Godunov scheme type 
numerical method was developed to solve the coupled set of 20 moment 
transport equations for a quasineutral single-ion plasma. The numerical 
method and the first results were presented at several national and interna- 
tional meetings and a paper describing the method has been published in 
the Journal of Computational Physics. To our knowledge this is the first 
numerical method which is capable of producing stable time-dependent solu- 
tions to the full 20 (or 16) moment set of transport equations, including the 
full heat flow equation. Previous attempts resulted in unstable (oscillating) 
solutions of the heat flow equations. Our group invested over two man- 
years into the development and implementation of the new method. The 
present model solves the 20 moment transport equations for an ion species 
and thermal electrons in a domain extending from a collision dominated to 
a collisionless region (200 km to 12,000 km). This model has been applied 



to study 0+ acceleration due to Joule heating in the lower ionosphere. 

Analytic investigations of transport phenomena. Another basic element of 
our proposed research program was the investigation of basic transport prop- 
erties of various plasma components in different regions of the solar system. 
These studies led to several interesting results. A full analytic solution (den- 
sity, flow velocity, parallel and perpendicular pressure) was found to the 
double adiabatic equations in a rotating planetary magnetosphere. Approx- 
imate analytic solutions were found for the heat flow equations in the topside 
ionosphere. A second-order quasilinear theory was developed to model the 
effects of wave-particle interactions in geoplasmas. We also developed a new 
2D model of flux tube convection and used it to explain unusual phenomena 
observed during a huge geomagnetic storm. 

Modeling of polar wind and plasmaspheric flows. With the help of our SPTP 
grant we developed a time-dependent, multispecies hydrodynamic model of 
interhemispheric plasma flow for tilted geomagnetic field lines. The model 
flux tube connects two reservoirs, each at an altitude of 200 km; these reser- 
voirs represent photochemically controlled regions of the ionosphere. The 
model successfully describes high speed plasma flows between two photo- 
chemically controlled conjugate ionospheres. We also developed new models 
of superthermal electron transport in the ionosphere and plasmasphere. 

Studies of the cometary pick-up process. The pickup of freshly ionized par- 
ticles in cometary environments turned out to be one of the most fascinating 
problems of theoretical space plasma physics. The lessons learned at comets 
can also be applied extensively in a broad range of other areas (such as 
interstellar pickup, solar wind termination, shuttle outgassing, active exper- 
iments, barium releases, planetary magnetospheres). In the last couple of 
years we addressed several aspects of the cometary pickup problem: wave 
generation by pickup ions, kinetic description of particle pitch-angle scat- 
tering, and first and second order Fermi acceleration of pickup ions. A new 
transport equation was also derived which preserves causality in the trans- 
port of pickup ions. 
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